Desorption of alkali atoms from 4 He nanodroplets 
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The dynamics following the photoexcitation of Na and Li atoms located on the surface of helium 
nanodroplets has been investigated in a joint experimental and theoretical study. Photoelectron 
spectroscopy has revealed that excitation of the alkali atoms via the (n + l)s <— ns transition leads 
to the desorption of these atoms. The mean kinetic energy of the desorbed atoms, as determined 
by ion imaging, shows a linear dependence on excitation frequency. These experimental findings 
are analyzed within a three-dimensional, time-dependent density functional approach for the he- 
lium droplet combined with a Bohmian dynamics description of the desorbing atom. This hybrid 
method reproduces well the key experimental observables. The dependence of the observables on 
the impurity mass is discussed by comparing the results obtained for the 6 Li and 7 Li isotopes. The 
calculations show that the desorption of the excited alkali atom is accompanied by the creation of 
highly non-linear density waves in the helium droplet that propagate at supersonic velocities. 

PACS numbers: 36.40.-c, 32.30.Jc, 78.40.-q, 67.40.Yv 
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I. INTRODUCTION 

During the last decade, many properties of helium 
droplets have been disclosed, especially in relation to 
their use as a gentle matrix for spectroscopic experi- 
ments. In particular, the dynamical aspects associated 
with the vibrational and rotational degrees of freedom 
of embedded chromophores, as probed by infrared spec- 
troscopy, are nowadays well established. We refer the 
interested reader to a series of review papers devoted to 
this subject^ - — At variance, the effect of the strong per- 
turbations induced by electronic excitation or ionization 
of impurities in helium droplets is much less understood. 
For example, there is still no consistent model that can 
fully account for the lineshapes and splittings observed in 
the electronic spectra of aromatic molecules^ Even less 
is known about the ensuing dynamics, like the rearrange- 
ment of the helium or the evolution of photoelectrons^iS. 
These aspects are however expected to be relevant for the 
study of chemical reactions in this soft ultra-cold envi- 
ronment, since chemical reactions involve essentially the 
rearrangement of the electronic structure of reactants. 
In this respect, understanding the dynamical evolution 
of an electronically excited impurity in helium droplets 
might be considered a first step towards a better under- 
standing of chemical reactions in this unique quantum 
environment. 

The aim of the present work is to gain insight into the 
dynamics initiated by the excitation of alkali atoms resid- 
ing on the surface of helium droplets. The np <— ns tran- 
sitions of alkali-doped helium nanodroplets have been the 
subject of a series of experimental and theoretical stud- 
ies, see e.g. Refs. IllTllol and references therein. Only 
recently, studies involving transitions to higher excited 



states have been reported These studies reveal that 
excitation of the alkali atoms in almost all cases leads 
to the desorption of the excited atoms from the surface 
of the helium droplet on a picosecond timescale i 19 ^ 20 In 
the present work we investigate, both experimentally and 
theoretically, the desorption process of Na and Li atoms 
following excitation to the 4s and 3s state, respectively. 
The use of two different atomic systems allows us to iden- 
tify general processes, while the use of lithium having 
two light isotopes, 6 Li and 7 Li, allows us to quantify the 
mass effect on these processes. In the experiments we 
have used ion and electron imaging techniques to deter- 
mine the state distributions of the desorbed atoms and 
the velocity distributions of these atoms as function of 
excitation frequency. The functional correlation between 
excitation energy and final kinetic energy of the impuri- 
ties has also been determined from first principles using 
Fermi's Golden Rule. The dynamics of the desorption 
process has been modeled using a newly developed the- 
oretical hybrid approach based on time-dependent DFT 
for the helium density and quantum trajectories for the 
impurity wave function. 

The hydrodynamic formulation of quantum mechan- 
ics -first developed by Bohm-^i has drawn the atten- 
tion of theoretical quantum chemists and physicists be- 
cause it constitutes an alternative representation to the 
time-dependent Schrodinger equation that allows to over- 
come some of the computational problems inherent to 
the conventional quantum mechanical approach^ The 
Bohmian formulation bears a large flexibility, and a va- 
riety of quantum trajectory methods are found in the 
literature, each of them adapted to the nature of the dy- 
namical problem under investigation^"^ We will show 
how the Bohmian method also lends itself to the study of 
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the complex dynamical processes addressed in the present 
work. 

This paper is organized as follows. In Sec. II we de- 
scribe the experimental setup. The quantum-trajectory 
model developed here to analyze the experimental results 
is presented in Sec. III. In Sec. IV the experimental re- 
sults are presented and compared with the outcome of the 
calculations. Further remarks on the results are made in 
Sec. V, and we conclude with a summary, Sec VI. Finally, 
a brief theoretical discussion on the angular dependence 
of the velocity distributions is presented in Appendix A, 
and the expressions used to simulate the experimental 
observables are collected in Appendix B. 



II. EXPERIMENTAL 

The experimental setup has been described in detail 
before i 27 i 28 In brief, helium droplets consisting on av- 
erage of several thousands of atoms are formed by ex- 
panding high-purity helium gas at a pressure of 30 bar 
into vacuum through a 5 /xm orifice cooled to cryogenic 
temperatures. The size distribution of these droplets 
can be systematically varied by changing the source 
temperature^ The helium droplets pick up alkali atoms 
as they traverse an oven in which either sodium or lithium 
metal is evaporated. The temperature of the oven is ad- 
justed to ensure that the droplets on average pick up 
less than one alkali atom. Via a differential pumping 
stage the doped droplets enter a velocity map imaging 
spectrometer. At the center of this setup, the alkali- 
doped droplets are excited by crossing the droplet beam 
perpendicularly with the frequency-doubled output of 
a Nd:YAG pumped dye laser (PrecisionScan SL, Sirah 
Laser- und Plasmatechnik GmbH). The laser system is 
operated at a repetition frequency of 20 Hz and provides 
radiation with a linewidth of less than 0.1 cm -1 , an en- 
ergy of 5 mJ/pulse and a pulse duration of 11±1 ns. The 
laser beam is slightly focused to yield an estimated spot 
size of 0.37 mm 2 at the excitation region. Following exci- 
tation, the alkali atoms are ionized by the absorption of 
an additional photon from the same laser pulse. The ions, 
or alternatively the photoelectrons, are accelerated by the 
applied electric fields and projected onto a position sen- 
sitive detector consisting of a pair of microchannel plates 
and a phosphor screen. The light emitted by the phos- 
phor screen is imaged onto a high-resolution CCD camera 
(A202k, Basler) that is read out every laser shot. The in- 
dividual images are analyzed online and the centroids of 
the impacts are determined. The kinetic energy distri- 
butions are determined by performing an inverse Abel 
transform on the image constructed from the accumu- 
lated centroids. Spectra are recorded using this setup by 
monitoring the number of ion impacts on the detector 
as function of laser frequency. Both the ion images and 
the spectra can be recorded at a specific mass by gating 
the front of the detector at the arrival time of the ions of 
interest. 



III. THEORETICAL APPROACH 

The theoretical model is described here in detail for the 
4s <— 3s excitation of Na impurities on 4 He droplets, but 
the formalism applies equally well to the corresponding 
3s <— 2s transition in lithium. We would like to point 
out that, while the (n + l)s <— ns transitions are dipole 
forbidden in the gas phase, they are allowed when the 
atoms reside on the surface of helium droplets due to the 
reduced symmetry of the systemJi 

Our starting point is the Fermi Golden Rule (GR) for 
optical transitions derived from perturbation theory^ 
It yields the transition probability per unit time 
from an initial state \i) to a final state |/) due to 
the interaction with a perturbative electromagnetic field 
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where T-L is the unperturbed hamiltonian that describes 
the system. For our purpose here, the initial state is the 
ground state (gs) of the complex consisting of a superfluid 
4 He droplet and a Na atom at zero temperature. We 
label it as \i) = \^^ c , ip^ s a , <pl s ), where ^Hc is the nuclear 
many-body wave function of the helium cluster, i/j^a the 
nuclear wave function of the sodium atom, and ip e the 
electronic wave function of the complex. For convenience, 
we make explicit that the Na valence electron is in the 
nominal 3s state, as it is this electron that interacts with 
the electromagnetic field. We are interested in those final 
states I/} that are accessible by the optical transition to 
the nominal 4s state and we hence label them as |/} = 

To evaluate the integrand of the GR we make the fol- 
lowing approximations: 

i) the Born-Oppenheimer (BO) approximation 30 is 
used to factorize the electronic and nuclear wave- 
functions. 

ii) Density functional theory (DFT) 4 ^ is used to de- 
scribe the droplet-impurity complex, factorizing 
the Na and He nuclear wavefunctions, both evolv- 
ing according to self-consistent mean-field hamilto- 
nians. 

iii) the Franck-Condon (FC) principle 3 - - is invoked, so 
that it can be assumed that the atomic nuclei do 
not change their positions or momenta during the 
electronic transition. 



A. The electronic contribution 

The BO approximation allows one to factorize 
the electronic contribution in the GR expression as 
|(</? 4s |V/|<,9 3s )| 2 , and to use effective pair-potentials for 
the nuclear hamiltonians. Due to the dipolar nature of 
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the transition, the angular dependence of Kv^l^lVe")! 2 
can be written as P(9) = ^[1 + f3 P 2 (cos 9)} where 
P%(x) is the second Legendre polynomial, 9 is the angle 
between the direction of the polarization vector of the 
laser light and the final velocity of the adatom, and /3 is 
the so-called anisotropy parameter. 

Since the projection of the orbital angular momen- 
tum onto the symmetry axis of the system — defined 
by the sodium atom and the center-of-mass of the he- 
lium droplet — does not change during the electronic 
transition, the value of /3 can be inferred by considering 
the symmetry of the valence electron wavefunctions, ipf* 
and iff 1 . Under the assumption that both these nom- 
inally spherical wavefunctions exhibit a dipolar defor- 
mation and that this deformation is larger for the final 
wavefunction ip^ s as a result of the stronger interaction 
with the helium cluster,— one finds that the value of the 
anisotropy parameter is limited to 1.5 < j3 < 2, see Ap- 
pendix A. 



B. The nuclear contribution within DFT 

For the nuclear contribution, we start from the DFT 
description of the system in its ground state as described 
e.g. in Refs.[3and[3ll. We want to recall here that this is 
a zero temperature phenomenological description of su- 
perfiuid liquid 4 He that reproduces its thermodynamic 
properties and elementary excitation spectrum. The BO 
factorization of the electronic wavefunction allows one to 
represent the interaction between the helium moiety and 
the impurity by an effective He-Na interaction which is 
based on the Vx s (rH C -Na) pair-potential, rjje-Na being 
the interatomic distanced From the minimization of the 
energy density functional E[^/n c , V'Na], we obtain the ef- 
fective hamiltonians corresponding to both the helium 
particle density /9Ho( r ) = I^HcMI 2 , and the Na wave- 
function V'Na(r): 

H^ He (r) = |-^ + f/ [p Hc ] + ^H 3 e S (r)}*Hc(r) 

= Ai*Ho(r) 
h 2 \7 2 



2m Na 
essV'Na(r) 



+ V&(t) ^Na(r) 
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Here, f/[pHe] is the helium DFT self-consistent potential, 
/j, the helium chemical potential and e gs the ground state 
eigenenergy of the sodium atom in the mean field cre- 
ated by the helium droplet. The mean field interaction 
potentials are obtained by convolution 
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V*l(r) = J dr'|* H c(r')| 2 V^(|r'-r|) 



(3) 



These density-dependent, mean field hamiltonians are 
employed to find the ground state of the droplet-impurity 



complex — \&jj e and tf^ — used in the evaluation of the 
GR. Within the FC approximation the helium density is 
frozen, which allows us to identify with VE"^* and to 
write the total transition probability as 
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with 
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(5) 

and H^f a is the effective hamiltonian 
for the excited impurity, calculated by evaluating Eq. ([3]) 
using the ground state helium density and the the excited 
state pair-potential V" 4s (rHc-Na) of Ref. l35l. Introducing 
the identity ^ n |n)(n| = 1, where {|n)} is a complete set 
of Hyf a eigenstates with eigenenergies fkj n , we obtain 
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It is worth noting that summing Eq. (|6]) over the final 
states / one obtains the Franck-Condon factors for the 
Na absorption spectrum^Sr— 
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As previously indicated, the final states ^ a are those 
accessible by the optical transition to the nominal 4s 
state. By using the same |n) basis for the final states 
we recover from Eq. © the Franck-Condon factors. 
However, in the ion imaging experiments discussed be- 
low what are probed are adatoms after their desorption 
from the helium cluster, which are characterized by a lin- 
ear momentum k. For this reason, we are interested in 
those I/} states that, after a time evolve to asymptot- 
ically free states with well-defined momentum k, namely 
I^Na) = ^Na(0,too)|k), where J7 Na (0, is the quantum 
time-evolution operator. 

In this evolution we assume that i) no Na-He exciplexes 
are formed during the desorption process, and ii) the he- 
lium does not induce relaxation of the excited adatoms. 
These assumptions, which are justified a posteriori by the 
experiments discussed below, imply that all the states ac- 
cessible in the course of the excitation evolve to free 4s 
states of Na. This ensures that the free 4s states of Na, 
|k), represent a complete basis. Introducing these states 
into Eq. © we obtain 
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Integrating Eq. ([5]) over to one obtains the Na probability 
density in momentum space after the desorption process 



J(k) = J dw/(w,k) 
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which is a measured observable in our experiments. It 
is worth noting that the evolution defined by f (0, too) 
describes the processes that follow the optical excitation. 
Hence, at this point the FC approximation is no longer 
applicable and the helium density is allowed to evolve 
dynamically. 



C. ID exploratory calculations 

To obtain the dynamical evolution, we have to solve 
the coupled 3D time-dependent system 

^JU Ho (i,r) = fl&(t)*ito(i,r) 

*ftJ^N.(t,r) = fl&(t)^Mt,r) (10) 

from t = to too, using as initial condition the 3s ground 
state for both the helium and the Na nuclear wave func- 
tions. 

Before attempting the solution of Eqs. (|10p. we have 
carried out an exploratory ID evolution. Firstly, we have 
solved Eqs. © by means of imaginary time methods -see 
e.g. Refs. I36U37I for details- obtaining the 3D structure of 
the ground state of a Na@ 4 Heiooo droplet (shown in the 
top left panel of Fig. H~2j) . Secondly, keeping the helium 
density frozen, we have let the Na wave function evolve 
in the resulting mean field potential along the symme- 
try axis, z, defined by the center-of-mass of the helium 
droplet and the sodium atom. Since the electronic exci- 
tation is assumed to be instantaneous in the FC approxi- 
mation, the calculation starts from the Na wave function 
corresponding to the ground state of the V^l potential, 
which next evolves in the potential shown in Fig. [1] 
We solve the Schrodinger equation for Na in a regular 
mesh using a predictor-corrector algorithm, 39 whose first 
time steps are provided by the outcome of a fourth-order 
Runge-Kutta algorithm. We have found that the Na 
atom reaches a quasi-free motion regime after about 1 
ps. Since the helium density is not allowed to evolve, the 
very large mass of the helium droplet causes that all the 
potential energy deposited into the system during the ex- 
citation is converted into kinetic energy of the Na atom. 
This is in clear disagreement with the experimental re- 
sults discussed below that show that a significant amount 
of energy is transferred to the helium droplet. The ID 
calculation is nonetheless useful as it helps to identify 
limitations of the approach and it enables us to deter- 
mine time and space scales for the full 3D simulations. 

Due to the highly repulsive nature of the VJja poten- 
tial, the Na atom leaves the droplet very quickly and 
attains a high mean asymptotic velocity of ~ 650 m/s. 
As a result of the high velocity, fast oscillations appear 
in the Na wave function as shown in Fig. [5J In order to 
reproduce these high frequency oscillations the use of a 
very fine spatial mesh is mandatory. This makes a full 



3D evolution computationally unaffordable. Fortunately, 
both the probability density and the velocity field of Na 
are smooth functions which makes it possible to describe 
the full 3D dynamics of the Na atom with quantum tra- 
jectories, an approach that uses positions and velocities 
instead of complex wave functions. This allows us to use 
a hybrid calculation scheme in which we compute the he- 
lium wave function in a mesh using standard methods for 
partial differential equations, 39 while the evolution of the 
Na atom is solved using Bohmian dynamic s 21 ' 22 as indi- 
cated below. We mention the existence of other hybrid 
simulations for the description of the dynamics of doped 
helium droplets 42r— 



D. Bohmian trajectories for the impurity dynamics 

The equation of motion for the Na trajectories R(t, r) 
is derived as follows i 22 i 24 ~ — We begin with writing 
the Na wave function in its polar form ^N a (i, r) = 
^K{t,r)e iS( - t ' r)/n , where 7£(t,r) is the probability den- 
sity and S(t,r) is the phase in units of h. S(t,r) is also 
known as the velocity potential since the velocity field 
is defined as v(i, r) = V S(t,r)/m^ a . Splitting the Na 
time-dependent Schrodinger equation into its real and 
imaginary parts, we obtain a continuity equation for the 
probability density coupled to a Hamilton- Jacobi (HJ) 
equation for the phase 



8U{t, r) 



= V-j(t,r) 



= \m^\w{t, r)| 2 + Q(t, r) + V&(t, r) , 

(11) 

where j(t, r) = lZ(t, r)v(t, r) is the current density and 



Q(t,r) 



2rriNa y/K(t, r) 



(12) 



is the so-called quantum potential. 

Equations (fTTT) have been solved as follows. Writing 
the density and the current density at time t as an his- 
togram of M test particles with trajectories {Rj(i)} 
where Rj(i) = R(t, r,;) and R;(0) = r$, we get 



M 



1 M 

K(t,r) = lim ^[r-R^i) 

M— *-oo M z — ' 
M 



i=l 



j(t,r) = lim jTvtRiWMr-RiW] 



(13) 



The continuity equation is automatically fulfilled if 
Ri(£) = v[Rj(i)], i.e., if the change in time of the posi- 
tion of the test particle is just the velocity field evaluated 
at the position of that test particle. The equation of mo- 
tion obeyed by the velocity field -and thus the equation 
for the trajectories- is obtained by taking the gradient 
of the HJ equation and rewriting it in the Lagrangian 
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reference frame (d/dt = d/dt + v • V). One then obtains 
the quantum Newton equation 

m N aRi(*) = - V [Q(t, r) + V*t(t, r)] | r=Ri(t) . (14) 

In this way, both helium wave function and Na trajecto- 
ries are computed consistently at each time step of 10 -4 
ps using a fourth-order Runge-Kutta algorithm. The 
quantum potential Q(t,r) is computed using the his- 
togram of the test particles as probability density in a 
regular mesh with a spatial resolution of 0.35 A, and 
13-points formulas for the derivatives involved. Simul- 
taneously, the helium wave function is evolved with a 
fourth-order Runge-Kutta algorithm, using as initial con- 
dition the 3s ground state helium density and a set of 
M = 200 000 positions randomly generated from the 3s 
ground state Na probability density. The system is solved 
up to too = 5 ps (too = 3 ps in the case of Li) , when the 
test particles are far enough from the droplet to follow 
quasi-free trajectories at constant velocity. Note that this 
time is five times longer than that of the ID exploratory 
calculation with the frozen helium density. 

E. Practical evaluation of the Fermi Golden Rule 



the GR is semiclassical, the trajectories are quantum- 
mechanically determined. 

Incorporating shape fluctuations of the helium droplet 
around the impurity has proved to be crucial for achiev- 
ing a quantitative description of processes involving im- 
purities in helium, as they substantially contribute to 
the broadening of observable a 14 ' 45 i 46 . For this reason, we 
have included the effect of fluctuations in the evaluation 
of Eq. (|16p a posteriori, i.e., after the time evolution, us- 
ing the DF-sampling method, as shown in the Appendix 
B, Eq. (IB5p . Inclusion of the density fluctuations mainly 
affects the width of physical observables such as the ex- 
citation spectrum and the velocity and kinetic energy 
distributions of the desorbed atoms. 

We want to stress that the transition probability u>i->/ 
contains all the physical information that can be exper- 
imentally determined. In point of fact, i) by integrat- 
ing over k we obtain the excitation spectrum, as seen in 
Eq. (JT)); ii) by fixing the excitation energy uj, we obtain 
the distribution over momentum k, corresponding to the 
velocity distribution of the desorbed atom; and iii) by the 
change of variable Ek = fi 2 k 2 /2m^ & the kinetic energy 
distribution is obtained. The explicit expressions used to 
compute these quantities are reported in Appendix B. 



To complete the evaluation of the GR expression Eq. 
©, we have resorted to a semiclassical approximation for 
the i?Na hamiltonian for obtaining the dipole absorption 
spectrum of atomic impurities^ the kinetic term is ne- 
glected and the hamiltonian is replaced by the potential 
energy, H$f a (Q) — > V^l(0,r). Thus, its eigenstates are 
those of the position operator, \n) — > |r), and its eigen- 
values those of the potential energy surface evaluated at 
r, thus huj n — > 14 s (r) and the sum over states becomes 
an integral J2 n I allowing us to write 

J(w,k) = Jdr |<k|tf N B(t«»0)|r)| 2 

x |(r|C)| J 6 \f*> - O^s (r) - hu g ,)] . 

Within the quantum trajectory description for the Na 
atom, the probability |(k|E/Na(ioo, 0) | r ) | 2 is written as 
5 [k — K(too, r)], where K(t oc , r) is the momentum at too 
of the trajectory with initial position r at t = 0. Since 
Na moves at too as a free particle, we can safely con- 
sider KK(too,r) — mNaR-^oo, r). Using the definition in 
Eq. (fT3| for the Na probability density, Eq. (fl5|) is then 
computed as 

1 M 

I(u, k) = — ^ S^hk- m N aR. 4 (ioo) 

x8 {hu - [V&(0,Ri(p)) - hu ga ]} . 

The semiclassical approximation incorporates the trajec- 
tories in a natural way by correlating along each tra- 
jectory its initial potential energy with its asymptotic 
linear momentum. Note that though the evaluation of 



IV. RESULTS 
A. Excitation spectra 

To investigate the desorption dynamics of sodium and 
lithium atoms from the surface of helium nanodroplets, 
the alkali atoms have been excited to their nominal 4s 
and 3s state, respectively. The corresponding excitation 
spectra are shown in Figs. 3 and 4. The spectra of the 
two lithium isotopes, 6 Li and Li, have been recorded in- 
dividually by gating the detector at the appropriate ar- 
rival time of these ions. It should be noted that only bare 
sodium and lithium ions have been observed and that no 
complexes with helium were detected, as was explicitly 
checked by time-of-flight mass spectrometry. This im- 
plies that the excitation spectra reported in Figs. [3] and 
U] correspond to absorption spectra which allows for a di- 
rect comparison between the experimental and theoreti- 
cal spectra. The spectra are all characterized by a broad 
absorption band that is blue-shifted with respect to the 
transition in the free atom. The 4s <— 3s transition of 
sodium-doped droplets has been discussed in detail be- 
fore and the large blue-shift has been attributed to the 
repulsive character of the 4s effective potential,— see also 
Fig.m 

The spectrum of the 4s 4- 3s transition for the 
Na@ 4 Heiooo system has been calculated taking explic- 
itly into account dimple fluctuations using the atomic- 
like DFT sampling technique, Eq. ([B7) . The width of the 
experimental spectrum, 400 cm -1 , is well reproduced by 
the calculations which yield a value of 370 cm -1 . In con- 
trast, the spectral shift, which is droplet size dependent 
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and amounts to 570 cm -1 for droplets consisting on av- 
erage of 1700 helium atoms, is somewhat underestimated 
by the calculation that yields 350 cm -1 for N = 1000. 

We would like to mention that although the absorption 
spectra of Na*^ and other impuritie o 13 ! 36 calculated by 
semiclassical methods yield at times quantitative agree- 
ment with the experimental spectra, the DFT sampling 
method employed here constitutes a consistent frame- 
work that is able to reproduce the basic spectral features 
both in droplets and bulk helium i 14 i 45 ~ — It is evident 
that since the helium-helium correlations are described 
in a semiclassical way, the method still needs some im- 
provements. In spite of this limitation, it is the only 
workable method for incorporating density fluctuations 
within the DFT scheme, and it is for this reason that we 
use it in this investigation. 

Inspection of the spectra for the two lithium isotopes 
shown in Fig. |4] reveals that they are very similar but 
that the spectrum of 6 Li is slightly less blue-shifted than 
that of 7 Li. Solving the equivalent of Eqs. @ for both 
isotopes indicates that this difference is related to the 
ground state structure: as can be seen in Fig. [5] the 
lighter 6 Li isotope generates a less pronounced dimple 
structure and is more delocalized than the heavier 7 Li 
isotope. As a result, the lighter isotope probes the ex- 
cited state potential at larger distances corresponding 
to lower energies. It is interesting to note that similar 
zero point motion effects have been observed for Li on 
3 He and 4 He droplets.— The calculations accurately re- 
produce the width of the absorption line (540 cm -1 in 
the calculations vs. 530 cm -1 in the experiments) and 
the isotopic shift (relative difference of 13% in the cal- 
culations and 15% in the experiments). As discussed 
above, the calculations somewhat underestimate the ab- 
solute blue-shift of the spectrum. We would like to point 
out that since the helium-lithium interaction is isotope- 
independent, all the differences found in the ground state 
structure -and thus in the excitation spectra- arise from 
the kinetic energy term in Eqs. ([5]) for Li. It is thus 
a pure quantum effect (zero point motion) and conse- 
quently cannot be reproduced by the calculations if the 
impurity is included as an external field. 



B. Photoelectron spectra 

The desorption efficiency and the helium-induced re- 
laxation of excited alkali atoms have been investigated 
using photoelectron spectroscopy. Fig. |6] shows the 
photoelectron spectrum obtained following excitation of 
sodium-doped helium nanodroplets at a frequency of 
26316 cm -1 , corresponding to the maximum of the ab- 
sorption band. The spectrum is characterized by a strong 
peak at low photoelectron kinetic energy and a much 
weaker peak at higher energies. Based on the photon 
energy and the ionization potential of sodium^ the low 
energy peak can be readily assigned to gas phase sodium 
atoms in the 3p state. The peak at high energy corre- 



sponds to free sodium in the 4s state. Since the spectrum 
reveals no other peaks which could be assigned to excited 
sodium atoms attached to the helium droplets,— we con- 
clude that all excited atoms desorb from the droplets on 
the timescale of the laser pulse. 

The photoelectron spectrum gives the impression that 
the helium induces a strong relaxation of the excited 
sodium, since the majority of the sodium atoms are found 
to reside in the 3p state and not in the initially excited 
4s state. However, when interpreting the photoelectron 
spectrum one has to take into account that the 4s state 
has a short radiative lifetime^ and that the ionization 
cross sections depend strongly on the excited state of the 
sodium atom.] 49 ! 50 Unfortunately, it is not possible to ad- 
dress this issue by recording the corresponding photoelec- 
tron spectrum of free sodium as the 4s <— 3s transition 
is dipole forbidden for one-photon excitation in the free 
atom. The photoelectron spectrum therefore has been 
simulated using a rate equation approach. The model 
used is graphically depicted in Fig. [7] Following exci- 
tation of the surface bound sodium atoms to the nomi- 
nal 4s state they desorb from the helium droplets. The 
free atoms decay by spontaneous emission to lower lying 
states. This process is characterized by the Einstein co- 
efficients A^-vap and A^ p ^ s ^ During the laser pulse 
the excited atoms are ionized with an efficiency that is 
determined by the state specific ionization cross sections 
CT4 S and c ^ 49 ' 50 and the laser intensity. To model the 
photoelectron spectra we now assume that no relaxation 
is induced by the helium and that the atoms desorb in- 
stantaneously from the droplets. This latter assumption 
is justified by the fact that the radiative lifetime and the 
interaction time of the free atoms with the light pulse 
is much longer than the desorption time of the excited 
atoms, vide infra. It should be noted that with these 
assumptions the calculated photoelectron spectrum cor- 
responds to that of a free sodium atom resonantly ionized 
by a 1 + 1 photon excitation process via the 4s state. The 
set of coupled equations describing the ionization of the 
free sodium atoms has been solved numerically assum- 
ing that the 11 ns laser pulse can de approximated by 
a Gaussian distribution. The uncertainties of the vari- 
ous constants used, as well as laser power fluctuations 
have been included in the simulations to determine the 
uncertainty of the relative intensities in the calculated 
photoelectron spectrum. 

The theoretically and experimentally determined in- 
tensity ratios are reported in Table Q] As can be seen 
from the table, the large intensity of the 3p state is well 
reproduced by the calculations. This signifies that the 
3j? state is mainly populated by spontaneous emission of 
4s excited sodium atoms during the laser pulse. The 
small discrepancy between the theoretical and experi- 
mental results might be attributed to the large uncer- 
tainty associated with the <j^ s ionization cross section, 
which is close to its minimum at the excitation frequency 
used. Alternatively, it might indicate that some helium- 
induced relaxation takes places during the desorption 
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process. In view of the good agreement between experi- 
ment and model calculations, we assume that all sodium 
atoms leave the helium droplets in the initially excited 
4s state. Although no photoelectron spectra have been 
recorded for lithium, we presume that they behave sim- 
ilar to sodium, i.e. all the excited lithium atoms desorb 
from the droplets without undergoing helium-induced re- 
laxation. 



C. Velocity and kinetic energy distributions 

To obtain insight into the desorption dynamics the ve- 
locity distributions of the desorbed atoms have been de- 
termined. To this end ion images have been recorded at 
several frequencies within the absorption bands. Fig. [5] 
shows two ion images of sodium that have been recorded 
following excitation at the low and high frequency end 
of the 4s <— 3s excitation spectrum. Both images are 
characterized by a strong anisotropic angular and radial 
distribution, indicating that the desorbing atoms leave 
the droplets with a well-defined velocity distribution. By 
performing an inverse Abel transform to these images the 
speed distribution and the angular anisotropy parameter 
(3 have been determined. 

The resulting speed distributions, which are also shown 
in Fig. [5J are found to depend strongly on the excitation 
frequency. While excitation at 26100 cm -1 yields sodium 
atoms with a most probable speed of 440 m/s, excitation 
at 26600 cm" 1 yields significant faster sodium atoms hav- 
ing speeds of 695 m/s. Whereas the speed distributions 
depend on the excitation frequency, the angular distri- 
butions do not show such dependence. The anisotropy 
parameter (3 is found to be independent of the velocity 
of the desorbing sodium atoms and has a mean value 
of 1.81±0.06. Similar values for the anisotropy parame- 
ter are found for desorbed lithium atoms, see Table [TTJ 
It should be noted that the speed distributions and the 
values for the anisotropy parameters are found to be in- 
dependent of helium droplet size. This can be attributed 
to the local character of the interaction of the alkali atom 
with the helium droplet. 

The values of the anisotropy parameter found in the 
experiments are close to that expected for the parallel 
(ri + l)s <— ns transitions. As discussed before, the 
small deviation from this value might be due to helium- 
induced configuration mixing. Recently, Callegari and 
Ancilotto proposed a method to calculate the interaction 
potentials of alkali atoms on helium nanodroplets that ex- 
plicitly takes into account configuration mixing. 51 Using 
the expansion coefficients determined by this method^ 
evaluation of Eq. (|A3j) yields a value of (i = 1.99 
for the anisotropy parameter. Evidently, configuration 
mixing cannot account for the experimentally observed 
anisotropy. A reduction of the anisotropy parameter also 
can result if the rotational period of the helium droplets 
is comparable to the time scale of desorption^ Even 
though not much is known about the rotation of helium 



droplets, it is to be expected that the rotational period 
depends on the size of the droplets. Since the anisotropy 
parameter is identical for all three adatoms and does not 
depend on the velocity of the desorbed atoms nor on the 
droplet size, it is highly unlikely that droplet rotation 
is the cause for the reduced anisotropy parameter. It is 
more likely that density fluctuations of the helium in the 
proximity of the alkali atoms during the desorption pro- 
cess lead to off-axis motion of the alkali atom. As we 
show in Appendix A, the deformations induced by the 
helium fluctuations can indeed result in a decrease of the 
anisotropy value. An estimation of the order of magni- 
tude of that decrease demands an explicit calculation of 
the electronic configurations, which is beyond the scope 
of this work. 

In order to establish a relation between the kinetic 
energy of the desorbing atoms and the excitation fre- 
quency, the speed distributions have been transformed 
into kinetic energy distributions, see Fig. [5J and their 
average value and standard deviation have been deter- 
mined. Fig. [TUJ shows the average kinetic energy of des- 
orbed sodium atoms as function of excitation frequency, 
while Fig. [TT] shows those for the two lithium isotopes. 
For all three impurities the average kinetic energy shows 
a linear dependence on the excitation frequency. The 
data points therefore have been fitted to the following 
expression: 

(Ekin) = ri(Tiui ~ h^o) (17) 

where huo corresponds to the excitation energy yielding 
alkali atoms with zero kinetic energy, while r\ is a pro- 
portionality constant. The constants derived from the 
fits are reported in Table |TIJ The values of ujq are very 
close to the transition frequencies of the free atoms, i.e., 
25740 cm -1 and 27206 cm -1 for sodium and lithium, 
respectively—^, The good agreement between these val- 
ues indicates that the proportionality constant rj can be 
interpreted as the fraction of the available energy that 
is converted into kinetic energy of the desorbing alkali 
atom. The values of rj as determined in the experiment 
depend not only on the type of alkali atom but also on 
its mass, see Table [TTJ The lighter the atom, the larger 
the value of r\ and thus the larger the fraction of avail- 
able energy that is carried away by the departing atom. 
The standard deviation of the kinetic energy distribu- 
tion, indicated by the bars in Figs. [TU] and [TTJ also shows 
a linear dependence on excitation frequency and there- 
fore has been fitted to the same functional form as the 
average kinetic energy. In contrast to the average kinetic 
energy, the proportionality constant for the standard de- 
viation, ?7a, does not depend on the atom nor its mass, 
see Table HI 

More insight into the desorption process can be gained 
from the calculations. As can be seen in Fig.[JJ excitation 
of Na to the 4s state leads to a highly repulsive interac- 
tion between the sodium atom and the droplet which 
triggers the desorption of the impurity. Figure [T^] shows 
the evolution of the Na@ 4 Heiooo system after excitation 
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of the sodium atom. Inspection of this figure reveals 
that the desorption of the sodium atom is accompanied 
by the creation of density waves in the helium droplet. 
This indicates that a fraction of the energy deposited in 
the system by the optical excitation is transferred to the 
helium cluster. It can be seen from the figure that af- 
ter ^1.5 ps the Na atom has left the droplet surface and 
the helium starts filling the dimple. The calculations are 
stopped after 5 ps, when the sodium atom has reached 
an asymptotic mean velocity of 510 m/s. During the 5 
ps time propagation, the position of the center-of-mass of 
the droplet experiences a minute displacement of 0.1 A. 
Assuming a constant motion, the translational energy of 
the droplet is calculated to be only 0.7 cm -1 , indicating 
that most of the energy transferred to the helium cluster 
is converted into internal energy. 

In the calculations, the GR is evaluated according to 
Eq. (|B1[) using the initial position and final velocities 
of the Na test particles as input. This expression gives 
direct access to the final velocity distributions of the des- 
orbed alkali atoms. To obtain a correlation between the 
asymptotic kinetic energy distributions and the excita- 
tion energy, the GR is evaluated using Eq. (|B2|) . The 
results of the calculations are compared with experimen- 
tal kinetic energy distributions in Fig. UH It can be seen 
that as the excitation energy increases, the kinetic energy 
distributions shift to higher energies and broaden. While 
these trends are well reproduced by the calculations, the 
widths of the distributions are somewhat overestimated. 
This can be attributed to the semiclassical method used 
to simulate the density fluctuations, as discussed before. 

The full I(u>,Ek) distribution calculated for Na is 
shown in Fig. [TU] by a false color representation. The 
average kinetic energy of the desorbed sodium atoms de- 
rived from this distribution at selected excitation energies 
is also presented in this figure. The calculations show a 
strong correlation between the mean kinetic energy and 
the excitation energy, in agreement with experiments. A 
linear fit to the theoretical data points yields a slope 
0.58 which compares well to the experimental value of 
0.52. It is important to point out that it is essential 
to take into account the density fluctuations in order to 
achieve a good agreement with experiment. If the density 
fluctuations are not included in the calculations a slope 
of 0.76 is obtained. This leads to the conclusion, that 
the density fluctuations are responsible for the absorb- 
tion of a large part of the energy deposited in the sys- 
tem by the electronic excitation. 4 - Calculations for the 
smaller Na@ 4 He5oo complex yield results similar to those 
for Na@ Heiooo- This size independence is in agreement 
with experiment, and can be attributed to the local char- 
acter of the interaction of the alkali atom with the helium 
droplet. 

The calculations reveal that lithium behaves very sim- 
ilar to sodium, although there are some differences. Due 
to the lighter mass and a less pronounced dimple, Li 
atoms desorb faster from the droplets than Na atoms 
and reach the asymptotic free regime already after ~3 ps. 



Their asymptotic mean velocities are significantly higher 
than for sodium and are slightly different for the two Li 
isotopes, 6 Li being faster than 7 Li. The I(ui,Ek) dis- 
tributions for Li are shown in Fig. [TT] together with the 
experimental data. Also in this case the calculations dis- 
play a similar correlation between the mean kinetic en- 
ergy of the desorbed atom and the excitation energy. In 
agreement with experiment, see Table [TT1 slightly differ- 
ent slopes for the two isotopes are found, 0.80 for 6 Li 
and 0.76 for 7 Li, confirming the mass-dependent nature 
of the process. 

As a final comment, we want to point out that while 
the experimental mean kinetic energies vary essentially 
linearly with excitation energy, the theoretical values 
show some non-linearity, especially at higher excitation 
energies. We have found that the variation becomes more 
linear at these energies if we do not include helium den- 
sity fluctuations. This indicates that at least part of the 
non-linearity arises from the way fluctuations are han- 
dled. 



D. Helium Density Waves 

Whereas the experiments provide only information on 
the desorbing alkali atom, the calculations also provide 
insight into the dynamics of the helium droplet upon ex- 
citation of the adsorbed atom. As an example, in Fig. [T51 
the helium density profile of a Li@ 4 Heiooo droplet along 
the symmetry axis is represented as function of time. 
This figure, together with Fig.[T2"lfor Na, gives a pictorial 
yet quantitative representation of the dynamics triggered 
by the excitation of the alkali atom. In particular, they 
show the dramatic changes in the droplet density caused 
by the excitation and subsequent desorption of the im- 
purity. 

To establish the nature of the helium density waves 
created by the excitation, we concentrate on the simula- 
tions for Na, as they last for longer times. Fig. [T4l shows 
the evolution of the helium density profile of a Na-dopcd 
4 Heiooo droplet during the first 5 ps. Initially, the droplet 
extends along the z (symmetry) axis from about 21 to 
—25 A, and the Na atom is located in a dimple at the 
surface. Excitation of the sodium to the 4s state causes 
the dimple first to deepen due the highly repulsive nature 
of the He-Na 4s interaction. The associated compression 
of the helium last up to ~ 1 ps, as shown in the figure. 
Following this compression, the helium surface bounces 
back and the dimple starts being filled. The more distant 
part of the droplet at ~ —25 A remains unperturbed and 
at rest, indicating that energy deposited in the droplet 
leads almost exclusively to the excitation of its internal 
degrees of freedom and not to a translational motion of 
the droplet as a whole. 

Figures [T2lfT4l all reveal that excitation of the alkali 
atom launches highly non-linear density waves into the 
helium droplet. This becomes even clearer in Fig. [15] 
which shows the density variations along the symmetry 
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axis at different times. The waves created by the excita- 
tion propagate in the droplet at supersonic velocities. 55 
In the case of Na, see Figs. fT4land[T5| the hrst perturba- 
tion front, labelled as 1, moves at ^890 m/s and reaches 
the opposite 'edge' of the droplet in less than 5 ps. This 
perturbation generates carrier waves with a phase veloc- 
ity between 300 and 370 m/s, modulated by supersonic 
envelope fronts with growing intensity. The ones with 
highest intensity, labelled as 2, have a group velocity of 
^590 m/s. The origin of this modulation can be traced 
back to the original structure of the droplet, being an 
'echo' of the solvation shells around the impurity in its 
ground state. Next, a high intensity wave appears trav- 
elling at ^370 m/s (labelled as 3), which generates sec- 
ondary waves propagating backwards. A closer analysis 
of this wave has allowed us to identify it with a solitary 
wave or soliton. Indeed, it can be seen from Figs. [T4l and 
[15] that it always corresponds to the maximum intensity 
of the travelling wave. 

We have found similar density waves in the case of Li, 
for which we recall that the evolution is stopped after 3 
ps. This time is long enough to establish that the first 
perturbation front moves at about 750 m/s for both iso- 
topes. This speed is lower than in the case of sodium 
which might be attributed to the smaller amount of en- 
ergy that is deposited into the droplet. Waves 2 and 3 de- 
velop as well, although they are not so clearly visible and 
their velocities cannot be determined with confidence. 

To our knowledge, the existence of these types of trav- 
elling waves has not been disclosed before. They bear 
some similarities with the waves produced in liquid he- 
lium by the de-excitation of electron bubbles ! 56 ' 57 These 
waves have been identified with shock waves in Ref. l56l . 
but with no clear justification for this statement. 



V. FURTHER REMARKS 

The most pertinent result of this study is undoubtedly 
the fact that the asymptotic kinetic energy of the desorb- 
ing alkali atoms scales linearly with the excitation energy. 
Both experiment and theory find that the amount of en- 
ergy carried away by the desorbing alkali atom depends 
not only on the atom but also on its mass. In particu- 
lar, the larger the mass of the impurity, the larger the 
amount of energy deposited in the helium droplet. 

Even though the calculations reproduce the experi- 
mental observations fairly well, they do not provide the 
physical insight required to identify the precise mecha- 
nism giving rise to the particular partitioning of the avail- 
able energy between the desorbing alkali atom and the 
helium droplet. Rather than a limitation of the method 
used, this is a signature of the actual complexity of the 
desorption process. 

Nonetheless, one would like to have some insight into 
the physical processes leading to the particular energy 
partitioning. In view of the analogy, one might consider 
describing the desorption of an excited alkali atom from 



the surface of a helium droplet as a photodissociation re- 
action. Various models have been put forward to describe 
the photodissociation of polyatomic molecules, each of 
them focussing on specific aspects of the process^ In 
case of a direct dissociation via excitation to a repulsive 
state, the process is best described by the simple im- 
pulsive model introduced by Busch and Wilson in their 
seminal work on the photodynamics of N02^ The ba- 
sic idea of the impulsive model is that the force needed 
to break apart the molecule is solely directed along the 
axis of the dissociating bond. Under the assumption that 
the bond breaks promptly, the available energy is ini- 
tially partitioned, by conservation of linear momentum, 
between the kinetic energies of the two atoms forming the 
bond. The kinetic energy of the atoms is subsequently 
partitioned between the translation, rotational and vi- 
brational degrees of freedom of each fragment. This im- 
pulsive model might be well suited to describe the des- 
orption of an alkali atom from the surface of a helium 
droplet given that excitation of the alkali atom leads to 
a highly repulsive interaction between the excited atom 
and the helium droplet, see Fig. [TJ and a fast desorp- 
tion of the alkali atom. The main difficulty in applying 
the impulsive model to the helium droplet system lies in 
that the alkali atom is not bound to a single helium atom 
but interacts with all the atoms making up the droplet. 
However, due to its short range character, the repulsive 
interaction is by far the strongest with the nearby he- 
lium atoms located at the surface of the dimple^ 5 - This 
is also born out by the calculations which reveal that 
predominantly these helium atoms are displaced imme- 
diately after excitation of the alkali atom, see Fig. [12] 
One might therefore consider reducing the problem to a 
pseudo-polyatomic model in which the alkali atom having 
a mass of niAk is considered to be bound to the helium 
droplet via a single helium moiety with an effective mass 
of ra e ff that is to be determined. Assuming that the im- 
pulsive model can be applied to this strongly simplified 
description of the system, the kinetic energy of the des- 
orbed alkali atom, £^i„(Ak) is related to the available 
energy Ahuj according to: 

Skia(Ak) = — Ahuj . (18) 

TO off + TOAk 

The model thus yields a linear dependence of the ki- 
netic energy on excitation frequency, as is observed in 
the present study. The slope is directly related to the ef- 
fective mass of the helium moiety with which the excited 
alkali atom interacts. Since the alkali atom interacts with 
several helium atoms located at the surface of the dimple, 
we expect m e g to be larger than the mass of a single he- 
lium atom. Although m e g cannot be calculated a priori, 
it can be determined by fitting the experimental slope 
of the kinetic energy on excitation frequency, see Table 
im One finds an effective mass of ~16 and ^24 amu for 
Li and Na, respectively. This would signify that the 3s 
electron of Li interacts on average with 4 helium atoms, 
while the 4s electron of sodium interacts with 6 helium 
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atoms. The difference in the number of interacting he- 
lium atoms is thought to reflect the difference in electron 
orbit radius and dimple structure. 

We would like to stress here that in spite of the fact 
that the impulsive model offers an explanation for some 
of the experimental and theoretical observations, it is 
only an approximative description that lacks any predic- 
tive power. The effective mass required for this model can 
only be determined a posteriori from the experimental or 
theoretical data. Due to the approximative nature of the 
model this effective mass should be interpreted with care. 
This is exemplified by the results for the two lithium iso- 
topes, where according to the impulsive model 6 Li inter- 
acts with more helium atoms than 7 Li, see Table HU This 
contradicts the conclusion based on the difference in the 
excitation spectrum of the two lithium isotopes and the 
calculations of the ground state structure, which indicate 
that the lighter 6 Li atom is located further away from the 
droplet surface and consequently interacts less with the 
helium than the 7 Li atom. 

As a final remark, the weak dependence of the stan- 
dard deviation of the energy distributions on the alkali 
atom or its mass points to the helium as the source of 
the broadening. As seen in the calculations, helium den- 
sity fluctuations cause some dispersion in the observables. 
Since these fluctuations are to a large extent independent 
of the impurity attached to the droplet, it is expected 
that so is the dispersion of the observables. 

VI. SUMMARY 

We have carried out a combined experimental and the- 
oretical investigation of the desorption of Na and Li alkali 
atoms from the surface of helium droplets following ex- 
citation via the (n + l)s 4— ns transition. These systems 
are well suited to gain insight into the dynamics of this 
complex phenomena, since neither exciplex formation nor 
helium-induced relaxation of the impurity obscure the 
analysis of the experimental findings or its theoretical 
interpretation Additionally, the use of Li allows to ad- 
dress the mass effect in the desorption process by making 
a direct comparison between the results for 6 Li and 7 Li. 

The analysis of the experimental results has been car- 
ried out within a full dynamical, three dimensional ap- 
proach that combines a time-dependent DFT description 
of the droplet with a Bohmian description of the impu- 



rity. To the best of our knowledge, this is the first time 
that such a theoretical framework has been developed 
and applied to the desorption of impurities from helium 
droplets. 

The experiments reveal that the (n + l)s <— ns tran- 
sitions of Li and Na atoms located on helium droplets 
are significantly blue-shifted with respect to the corre- 
sponding gas phase transitions. They furthermore dis- 
close that excitation of the alkali atom leads in all cases 
to its desorption from the helium droplet and that the 
average kinetic energy of the desorbed atom depends lin- 
early on the excitation energy. These observations are 
all reproduced by the calculations, which allows us to 
have confidence in theoretical observations that cannot 
be experimentally verified. More specifically, the calcu- 
lations indicate that the energy deposited in the system 
by the excitation of the alkali atom leads to the creation 
of highly non-linear helium density waves that propagate 
through the helium droplet at supersonic velocities. One 
of such waves could be identified as a soliton. Based on 
the good agreement between experiment and theory it 
becomes also possible to identify the main physical in- 
gredients necessary for a quantitative description of the 
desorption process, namely: i) the zero-point motion of 
the impurities, ii) a full dynamical description of both the 
helium droplet and the impurity, and iii) the inclusion of 
helium density fluctuations. These concepts have been 
implemented in a numerical approach that is very robust 
and can be applied to the description of other dynami- 
cal processes involving atomic impurities propagating in 
helium dropletS i 27 i 28 ' 43 i 44 ' 60 Addressing exciplex forma- 
tion, as observed in several systemS ) 17 i 19 ' 41 i 42 ' 61 ~ — how- 
ever, still remains a challenge from a computational point 
of view that will likely require a test particle description 
of both the impurity and the helium droplet. 
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Appendix A 

We discuss in this Appendix the probability of an optical transition between mixed electronic orbitals. Defining 
the initial and final states in cylindrical symmetry as 

oo 
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(9) =^T / B e Y eo (9), 



1=0 



where Yi m are the spherical harmonics of order (£, m), the angular dependence of the dipolar transition probability is 
written as22 



where (£'010|^0) is a Clebsch-Gordan coefficient. Due to the properties of these coefficients we can write 

. \A e ^\ 2 e(3 + 21) + \A t + \\ 2 {2t 2 +1-1), 



M{l + t) - 3 



\Y m (9)\ 2 e^^^cosI 



where Pj(x) are the Legendre polynomials of order j. If the involved states are limited to £ < 1 

<Pi(6) = Y 00 +A 1 Y W (6) 
ip f (9) = Y 00 +B 1 Y 10 (6), 



we obtain for the transition probability 



^/(0) = ^ + ^cos 2 (0)xi-: 



-in 



1 +\A 1 \y(3\B 1 f) 



P 2 {cos9) = 1 + /3P 2 (cost 



(Al) 



(A2) 



(A3) 



It can be seen that the anisotropy can only take positive values. When \Bi\ 2 > |^4i| 2 , the anisotropy is constrained 
to 1.5 < f3 < 2, and if the initial state is spherical {A\ = 0) one finds the limiting value = 2. 

When non-axially symmetric deformations are included (to ^ 0) the transition probability for £ < 1 reads 



Pi_>/(0) = J2 \Bi m \ 2 \A e , m ,\ 2 \(e'm'lO\em)\ 2 \Y m (9)\ 2 cx 1 + f3 KP 2 (cos 9) 



(A4) 



with (3 = 2/(1 + |Aio| 2 /(3|Bio| 2 )) [same structure as in Eq. $K5§ ] and 



K 



(\A 10 \ 2 + 3|g 10 | 2 )(4|i? 10 | 2 - \A n \ 2 \B n \ 2 - 
2|5 1 o| 2 [2(|^ 10 p + 3|i3 1 o| 2 ) + 3(|A 11 | 2 |Bii| 2 + |A 1 _ 1 |2|S 1 _ 1 |2)]- 



(A5) 



Some algebra shows that K < 1, since (|Ai | 2 + 9|Bi | 2 )(|Aii| 2 |Bh| 2 + |^4i_i | 2 | 2 ) > 0. Consequently, these 
deformations will always decrease the value of the anisotropy parameter. 



Appendix B 



We show in this appendix how we have represented the delta functions involved in Eq. (Til)]) to compute the GR. 
Using the rectangular function rect(a;) = d{x + 1/2) — §{x — 1/2), where #(x) is the Step function, we write 



M 



Wi^f oc [1 + j3P 2 (cos 9)] — ^2 rect 



Ak 



1 

Ak' 



-rect 



nu>-[v£(o,B4(p))-nu gB ] 



Aw 



■L < (B1) 



where we have used for the number of test particles M = 200000 and the intervals take small values as Ak/m^ a ~ 5 
m/s and Aw ~ 3 K. The velocity distributions are simulated by choosing a value of the excitation energy w and 
evaluating Eq. (^Bl[) . while the final kinetic energy vs. excitation energy distribution is obtained after the change of 
variable E k = h k 2 /2m-^ a as 



I(oJ,E k ) 



2E k h 2 

M 



I \w,k 



2m Na E k 



— rect 



E k -TO N aR 2 (^oo)/2 

AE k 



AE k 



-rect 



ftw- [V^(0,Ri(0))-fiw ff ,] 



Aw 



1 
Aw 



(B2) 
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The excitation spectrum is computed as 



V ' M ^ 



rect 



Tkj- [V*l(0,Ri{0)) - hw g 
Aw 



1 

Aw' 



(B3) 



If density fluctuations are included using the stochastic method described in detail in Refs. HH the total distribution 
is generated by the contribution of N c configurations generated by sorting N random positions in the j-th configuration 
{ r n}n=i f° r the hard spheres that represent the N helium atoms, using the helium density divided by N as probability 
density distribution together with a hard-sphere repulsion between He atoms. The diameter of the sphere is of the 
order of the length h used to screen the Lennard- Jones potential and to compute the coarse-grained density£L This 
diameter is defined a a 14 ' 46 



1/3 



dl = h 



Po 



(B4) 



where po is the liquid saturation density value, and p is the coarse-grained density obtained by averaging the atomic 
density within a sphere of radius h£L The distributions are then obtained as 



JV C N 



M 



a jEE[ 1+ ^ ( cos *)] m ^ rect 

c j=l n=l i=l 



Ak 



Ak 



rect 



hu>- [^(1^-^(0)1)-^] 



Aa 
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Aw' 
(B5) 



- N a N ^ M 
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Ek — mNaR-f (*oo)/2 

~AE k 



AE k 



-rect 



hcj-[V^(\ri 



Ri(0)| 



Aw 



4~, (B6) 
Aw 



and 



J(w) = — V V — V rect 

7 — 1 n—l i=l 



Rw- [V^(|iJ -RiCO)!)-^,] 



1 

Aw' 



(B7) 



Note that with a number of configurations = 10000, the histograms are computed using a total of N c x M = 2 10 9 
contributions, large enough to reduce the statistical noise. ^ 
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FIG. 1: (Color online) Na initial probability density (thick 
solid line) used in the ID calculation and the 3s ground state 
(dashed line) and 4s excited state (solid line) mean-field po- 
tentials. The origin corresponds to the center-of-mass of the 
helium droplet. 




FIG. 2: (Color online) Real part of the Na wave function in 
the ID calculation at t = 0, 0.25, and 0.6 ps. Also shown are 
the probability density distributions (thick solid lines). 
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FIG. 3: (Color online) Theoretical and experimental exci- 
tation spectrum of the 4s <— 3s transition of Na on 4 Hejv 
droplets. The top-starred vertical line corresponds to the free 
atom transition. 48 
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FIG. 4: (Color online) Top panel: Experimental excitation 
spectrum of the 3s <— 2s transition for Li attached to helium 
droplets consisting on average of 6100 atoms. Bottom panel: 
Corresponding theoretical spectrum for Li@ 4 Heiooo- The top- 
starred vertical line corresponds to the free atom transition. 48 




FIG. 6: (Color online) Photoelectron spectrum recorded fol- 
lowing photoexcitation of sodium-doped helium droplets at 
26316 cm -1 . 




FIG. 7: Energy level diagram displaying the relation and exci- 
tation processes taking place after excitation of sodium-doped 
helium droplets via the 4s <— 3s transition, see text for details. 




21 






va4 ; 


1 1 ' 1 


26100 cm' 

26600 cm ' 


i . - . ,\ i 


I . 



200 400 600 800 1000 

Speed (m/s) 

FIG. 8: (Color online) Left: Ion images of sodium atoms des- 
orbed from the surface of helium droplets following excitation 
at two frequencies within the 4s -f- 3s resonance. The polar- 
ization of excitation laser is vertical with respect to the im- 
ages. Right: Speed distributions and anisotropy parameters 
derived from the ion images. 



22 



a 
3 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
Experiment 

/\ 

i i 
i ^ 

I 1 
I 1 


i i i 1 . i i i 1 i i . i 1 i . i 

26100 cm' 

26300 cm 1 

26600 cm 1 


,., .-i^'AV. , ,/x 





100 200 300 400 500 600 700 



3 



■e 

C3 



1 

Theory 


i 1 1 1 1 i 1 


i 


I I 

25900 cm"' 


/ X 
/ \ 






- 26100 cm 1 








26300 cm 1 


/ 

\ l 
i 

i \ 
i ^ 


\/ 




26600 cm 1 


i . 
i / 

/ 

i i 
/ 

/ 









100 200 300 400 500 600 700 
E kin ( cm ') 



FIG. 9: (Color online) Experimental (top panel) and theoret- 
ical (bottom panel) normalized kinetic energy distributions of 
desorbed Na atoms following excitation at different energies. 




FIG. 10: (Color online) I(u>,Ek) distribution for Na atoms. 
Squares: theoretical mean kinetic energy. Dashed line: lin- 
ear fit to the theoretical data up to an excitation energy of 
26400 cm^ 1 ^ Dots: experimental mean kinetic energy. Bars: 
experimental standard deviation of the kinetic energy distri- 
butions. Solid line: linear fit to the experimental data. 



24 




27 200 27400 27 600 27 800 28 000 28 200 28400 



Excitation Energy (cm 1 ) 
/ 




27 200 27400 27 600 27 800 28 000 28 200 28400 



Excitation Energy (cm 1 ) 

FIG. 11: (Color online) Same as Fig. [10] for Li. In this 
case, all the points shown have been included in the fit of the 
theoretical data. 
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FIG. 12: (Color online) Evolution of the Na@ 4 Heiooo complex 
after excitation. 
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FIG. 13: (Color online) Helium density profiles and Li prob- 
ability density distributions (Gaussian-like profiles) showing 
the dynamical desorption of Li isotopes along the symmetry 
axis. Dashed lines, 6 Li. Solid lines, 7 Li. 
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FIG. 14: (Color online) Evolution of the helium density profile 
of the Na@ 4 Heiooo system along the symmetry axis. Three 
supersonic fronts are identified and labeled by roman num- 
bers. Equidensity lines corresponding to 0.5 and 0.1 times the 
helium saturation density, 0.0218 A -3 , are shown in white. 
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FIG. 15: (Color online) Difference between the density profile 
and the initial density along the symmetry axis, Ap = p(t, z) — 
p(0, z), for t = 2.5 ps (dashed lines) and t = 5 ps (solid 
lines). The envelope wave of the modulated carrier waves 
is also shown for clarity (dotted lines). The three identified 
fronts are labelled as in Fig. 1141 
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State 


Experiment 


Model 


3 P 


0.96(1) 


0.89(4) 


4s 


0.04(1) 


0.11(4) 



TABLE I: Experimental and calculated relative intensities of 
the photoelectron peaks following excitation of sodium-doped 
helium droplets via the 4s <— 3s transition. 
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Atom 






Experiment 






Theory 




wo [cm ] 


V 




P 


m cf f [amu] 


ujq [cm j 


V 


6 Li 


27218(6) 


0.743(6) 


0.042(6) 


1.73(6) 


17.4 


27206 


0.802(8) 


7 Li 


27222(3) 


0.687(3) 


0.040(2) 


1.79(3) 


15.4 


27206 


0.756(8) 


23 Na 


25743(4) 


0.516(4) 


0.038(3) 


1.81(6) 


24.5 


25740 


0.583(9) 



TABLE II: Characteristics of the experimental and theoretical 
kinetic energy distributions of the desorbed alkali atoms, see 
text for details. 



